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On the possibility of empirically probing the Bohmian model in terms of the 
testability of quantum arrival/transit time distribution 

Dipankar Home* and Alok Kumar Pan^^ 
Department of Physics, CAPSS, Base Institute, Sector-V, Salt Lake, Calcutta 700091, India 

The present article focuses on studying the extent to which the nonuniqueness that is inherent in 
the standard quantum mechanical calculation of arrival/transit time distribution can be exploited 
to enable an empirical scrutiny of any causal trajectory model such as the Bohmian scheme. For this 
purpose, we consider the example of spin- 1/2 neutral particles corresponding to a wave packet which 

f^ , passes through a spin rotator(SR) that contains constant magnetic field confined within a region - 

in such a case, the transit time distribution can be measured in terms of the spin distribution of 
^ , particles emerging from the SR. In particular, we investigate the way one can compare the Bohmian 

predictions obtained for this example with that using one of the quantum approaches, say, the 
probability current density based scheme. Here the Bohmian calculational procedure involves a 

,S^ • couple of critical subtleties that lead to some specific directions for further studies. 

-I— ' 

C^ ■ I. INTRODUCTION AND MOTIVATION 

^. 

^- . . .2 

Born's interpretation of the squared modulus of a wave function (1-01 ) as the probability density of finding a 

'"^ particle within a specified region of space is a key ingredient of the standard framework of quantum mechanics, 

0^ thereby implying that the standard interpretation of quantum mechanics is inherently epistemological. On the other 

^ ... ... 2 

\^ hand, the possibility of an alternative interpretation of quantum mechanics by interpreting IM as the probability 

^D [2|developed the details of such an ontological model of quantum mechanics by using the notion of an observer- 

y—i , independent spacetime trajectory of an individual particle that is determined by its wave function through an equation 

^ ' of motion which is formulated in a way consistent with the Schroedinger time evolution. Bohm's model, thus, explicitly 

KS refuted the counterarguments (such as the ones put forward by Pauli[3[and von Neumann[4]) that claimed to have 

ruled out the formulation of such a model. Subsequently, much work has been done on various aspects of the Bohmian 

model[5]. That such a model is not unique has been elaborately discussed[6] while different versions of the ontological 
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model of quantum mechanics have been proposed 
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Although any such ontological model hinges on the notion of a definite spacetime track used to provide a description 
of the objective motion of a single particle, such trajectories are not directly measurable. Hence these trajectories 
have been essentially viewed as conceptual aids for understanding the various features of quantum mechanics. While 
recently a study has been reported which shows an application of such trajectories as computational aids for solving 
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the time-dependent Schroedinger equation[ll|here in this paper we explore the question as to whether an ontological 
model such as the Bohmian one is empirically falsifiable. 

The essence of the question as regards empirical equivalence between the standard and the Bohmian model can be 
seen as follows. Let us focus on the instant when an ensemble of particles begins to interact with a given potential. If at 
that initial instant, the position probability density of the ensemble is taken to be the same in both the standard and 
any ontological approach, then the time evolved position probability density of the final ensemble as calculated by the 
equation of motion specified by any ontological model is ensured to be the same as that obtained from the standard 
quantum mechanical formalism. The same equivalence holds good for any other dynamical observable represented by 
a Hermitian operator. 

The above is the main reason why the ontological models have been dismissed by many as merely of metaphysical 
interest, having a "superfluous ideological superstructure" (see, for example, the criticism of the Bohmian model by 
Heisenberg[l^- In this paper, we take a fresh look the question of empirical falsifiability of an ontological model like 
the Bohmian one. The motivation underlying the present study stems from the type of question as was posed by 
Gushing 1 13| : Could it be that the additional interpretive ingredients(such as the notion of particles having objectively 
defined trajectories) in the Bohmian model might enable testable predictions in a suitable example where the standard 
version of quantum mechanics has an intrinsic nonuniqueness allowing for a number of calculational approaches, but, 
prima facie, none of these can be preferred over the others using a rigorous justification based on the first principles? 
On this point one may recall that John Bell had once remarked that the Bohmian model of quantum mechanics is 
experimentally equivalent to the standard version "insofar as the latter is unambiguous" jl4| . 

It is from this perspective that the kind of nonuniqueness inherent in the standard quantum mechanical calculation 
of time distributions seems to be a pertinent tool for exploring the possibility of subjecting an ontological model such 
as the Bohmian one to an empirical scrutiny. 

Time plays a peculiar role within the formalism of quantum mechanics - it differs fundamentally from all other 
dynamical quantities like position or momentum since it appears in the Schrodinger equation as a parameter, not as 
an operator. If the wave function ip(x,t) is the solution of Schrodinger's equation, then 'ip(x,ti),ip{x,t2), il^{x,tj,). . 
. determine position probability distributions at the respective different instants ti, t2, t-^,. . . for a fixed region of 
space, say, between x and x + dx. Now, if we fix the positions atXi, X2, X3,. . ., the question arises to whether 
the quantities 'ip{Xi,t), 'ip{X2,t), ip^X^^t),. . . can specify the time probability distributions at respective various 
positions Xi, X2, ^3,. . .? However, one can easily see that, although J_ \%l:{x,ti)\'^dx — 1, one would have in 
general, J_ \ip{Xi,t)\'^dt ^ 1 . Hence, in order to quantum mechanically calculate the time probability distribution, 
unlike the position probability distribution, we do not readily have an unambiguously defined procedure. 

The fundamental difficulty in constructing a self-adjoint time operator within the formalism of quantum mechanics 
was first pointed out by Pauli|l5|. Another proof of the nonexistence of time operator, specifically for the time-of- 
arrival operator, was given by Allcock [l6|. Nevertheless, there were subsequent attempts to construct a suitable time 
operator. For instance. Grot et al. [l7| and Delgado et al. [ISjsuggested a time-of-arrival operator for a free particle, 
and showed how the time probability distribution can be calculated using it; interestingly, such an operator has an 



orthogonal basis of eigenstates, although the operator is, in general, not self-adjoint. 

In recent years there has been an upsurge of interest in analyzing the concept of arrival time distribution in 
quantum mecha nics, fo r useful reviews on this subject, see [19|,|20|. Here we note that a number of schemes have 



been formulated 
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-41| for calculating what has been called the arrival time distribution in quantum mechanics, for 
example, the probability current density approach, using the path integral approach, the consistent history scheme, 
and by using the Bohmian trajectory model in quantum mechanics. However, since there is an inherent ambiguity 
within the standard formalism of quantum mechanics as regards calculating such a probability distribution, it remains 
an open question as to what extent these different approaches can be empirically tested. 

There have been several specific toy models that have been to investigate the feasibility of how actually the mea- 
surement of a transit time distribution can be performed in a way consistent with the basic principles of quantum 
mechanics. The earliest proposal for a model quantum clock in order to measure the time of flight of quantum particles 



was suggested by Salecker and Wigner[42|, and later elaborated by Peres[43j.In effect, this model of quantum clock 
measures the change in the phase of a wave function over the duration to be measured. Such a model of quantum 
clock[43J can be used to calculate the expectation value for the transit time distribution of quantum particles passing 
through a given region of space. On the other hand, Azebel[44| has analyzed a process in which the thermal activation 
rate can serve as a clock. Applications of quantum clock models have also been studied for the motion of quantum 
particles in a uniform gravitational field by Davies |45fland others 35 1. 



Against the backdrop of such studies, in the present article we proceed as follows. Let us first consider the following 
simple experimental arrangement. A particle moves in one dimension along the x - axis and a detector is placed at 
the position x = X . Let T be the time at which the particle is detected, which we denote as the time of arrival of 
the particle at X. Can we predict T from the knowledge of the state of the particle at the prescribed initial instant? 

In classical mechanics, the time of arrival T at x — X oi an individual particle with the initial position xq and 
momentum po is T — t(t, xo,Po) which is fixed by the solution of the equation of motion of the particle concerned. But 
in quantum mechanics, this problem becomes nontrivial in terms of the probability distribution of times, say, H(r) 
over which the particles are registered at the detector location, say X , within the time interval, say, t and t + dt so 
that Jrp^ Il{T)dT is the probability that the particles are detected at a; = X between the instants Ti and T2. In other 
words, the relevant key question in quantum mechanics is how to calculate Il{T)dT at a specified position, given the 
wave function at the initial instant from which the propagation is considered. 

An effort along this direction was madef36] by considering the measurement of arrival time using the emission of a 
first photon from a two-level system moving into a laser-illuminated region. The probability for this emission of the 
first photon was calculated for this purpose by specifically using the quantum jump approach. Subsequently, further 
work|i37i]was done on this proposal using Kijowski's distribution[41]. 

In this work we address this question of empirical verifiability from a new perspective so that starting from any 
axiomatically defined transit/arrival time distribution one can directly relate it to the actually testable results. For 
this, we first need to evaluate the transit time distribution H(t) by using one of the suggested approaches. Using 
such a calculated this n(t), one can then derive a distribution of spin orientations n(0) along different directions 



for the spin-1/2 neutral particles emerging from a spin-rotator (SR) (which contains a constant magnetic field), (j) 
being the angle by which the spin orientation of a spin-1/2 neutral particle (say, a neutron) is rotated from its initial 
spin polarised direction. Note that this angle (j) is determined by the transit time (t)within the SR. Thus, in our 
scheme, the SR serves as a "quantum c/ocfc'Vhere the basic quantity is n(0)which determines the actually observable 
results which corresponds to the probability distribution of spin orientations along different directions for the spin-1/2 
neutral particles emerging from the SR. Such a calculated spin distribution function can be empirically tested by 
suitably using a Stern-Gerlach(SG) device, as explained later, thereby providing a test of n(<) based on which 11(0) is 
calculated. In this article we particularly focus on analysing the application of the Bohmian model in the context of 
such an example. 



3, 



In order to recapitulate the essential aspects of this setup |32| . we consider spin-1/2 neutral objects, specifically, the 
spin-polarised neutrons. The initial quantum state is given by '^ [x.t = Q) = tjj (x, t ~ Q)®x{t — ^) where x (t = 0) is 
the initial spin state taken to be the same for all neutrons whose spins are oriented along the x — axis. The spatial part 
is taken to be a propagating Gaussian wave packet (for simplicity, taken to be one dimensional) that moves along the 
-fir-axis, and passes through a region containing constant magnetic field(B), confined between x = Q and x = d{Fig. 
1). The enclosed magnetic field is along the -I- z— axis, and is switched on at the instant {taken to be t = 0) when the 
peak of the wave packet is at the entry point x — Q. 

In passing through this bounded region(henceforth designated as a spin rotator(SR)), the spin variable of a neutron 
interacts with the constant magnetic field, this interaction being mediated by the neutron's magnetic moment /i as 
the coupling constant. The angle (p by which the direction of the initial spin polarisation is rotated is determined 
by the time (t) over which this interaction takes place. Now, if the spatial wave packet is a superposition of energy 
eigenstates (plane waves), there will be a distribution(n(r)) of times r over which the interaction between neutron 
spins and the enclosed magnetic field takes place, corresponding to the relevant arrival/transit time distribution . 
Consequently, using the probability distribution function n(T), and an appropriate relation connecting the quantities 
(j) and T, one can calculate a distribution of spin orientations (11(0)) for the neutrons emerging from the SR. It is, 
therefore, evident that for calculating 11(0), the quantum mechanical evaluation of n(r) for the neutrons passing 
through the SR is a key issue. 

We will first elaborate on the relevant setup by discussing of how the determination of the quantity 11 (0) can 
actually be tested by using a SG device. Subsequently, we will discuss the key ingredients of the Bohmian procedure 
that can be adopted in the context of this setup for calculating 11(0) in terms of 11 (t). The subtleties involved in 
such a procedure will be critically analysed and directions for further studies will be indicated that are required to 
explore fully the potentiality of this example for subjecting the Bohmian model to empirical scrutiny. 

II. THE SETUP 

We consider an ensemble of spin 1/2 neutral particles, say, neutrons having magnetic moment /i. The spatial part 
of the total wave function is represented by a localised narrow Gaussian wave packet ^ {x,t — 0) (for simplicity, it is 
considered to be one dimensional) which is taken to be peaked at a; = at tha initial instant t = and moves with 




Figure 1: Spin-1/2 particles, say, neutrons with initial spin orientations polarised along the +x - axis and associated with a localized 
Gaussian wave packet (peaked at a: = 0, t = 0) pass through a spin-rotator (SR) containing a constant magnetic field B directed along 
the +z - axis. The particles emerging from the SR have a distribution of their spins oriented along different directions. Calculation of this 
distribution function is experimentally tested by measuring the spin observable along a direction n {9) in the xy-plane making an angle 6 
with the initial spin polarised along +x- axis. This is done by suitably orienting the direction n {0) of the inhomogeneous magnetic field 
in the Stern-Gerlach ISG — n\ device. 



the group velocity u. Thus the initial total wave function is given by ^ = -0 (x, t = 0) (g) x (^ = 0) where X (t = 0) is 
the initial spin state which is taken to be same for all members of the ensemble oriented along the x — axis. 

The SR used in our setup (Fig. 1) has within it a constant magnetic field B = Bz directed along the +z - axis, 
confined between x — and x = d. We consider a specific situation where the magnetic field is turned on while the 
peak of the initial wave packet reaches at the position x = at t — 0. 

Application of the Larmor precession of spin in a magnetic field has earlier been discussed, for example, in the 
context of the scattering of a plane wave from a potential barrier 46|, |47|. On the other hand, the scheme proposed 
here explores an application of Larmor precession such that one can empirically test any given quantum mechanical 
formulation for calculating the arrival/transit time distribution using Gaussian wave packet. 

Now, for testing the scheme for calculating the probability density function 11 (0) of spin orientations of particles 
emerging from the SR, we consider the measurement of a spin variable, say ct^, by a SG device (Fig.l) in which the 
inhomogeneous magnetic field is oriented along a direction n (6) in the xy-plane making an angle 9 with the initial 
spin-polarised direction (+x axis) of the particles. The initial x-polarised spin state can be written in terms of the 
z-bases | f ) z ^-'id | J, )^ as x(0) = 1/V2 (11)^ + 11 ) z)- While passing through the SR, the spin polarised state 
rotates only in the xy-plane by an angle (/) with respect to the initial spin orientation along x-axis. Such a rotated 



spin state in the xy-plane can be typically written as xi'P) — 1/V2 (| t 



rJ(t>\ 



■I) z)- If one applies the SG-magnetic 



field along the direction n (6), the relevant bases states corresponding to the spin operator ag are respectively | t)e = 
l/\/2 [\'\-)z+ e*^| i }z)j and | \)g = l/\/2 (| t )z - e'"] | )^). Then for such spin measurements, the probabilities 
of getting I t>0 and | i)g are p+(0) = |0(t | x(M^ - cos^B - 0)/2 and p_(0) = \e{i \ xW)\^ = sin^i^ - 0)/2 
respectively. 



III. THE TESTABLE PROBABILITIES 

Since we are considering an ensemble of spin-1/2 neutrons passing through the SR where initially all members 
of this ensemble have their spins polarised along, say, the +x axis, given the same initial spin state, the spins of 
individual members of the ensemble evolve over different times (characterised by 11 (i) , the distribution of transit 
times within the SR) of duration of the interaction with the constant magnetic field within the SR. 

Thus, for the spins of the particles emerging from the SR polarised along different directions (with the respective 
probabilities 11 (0) of making angles (p with the +x axis), the probabilities of finding the spin component along -\-6 
direction and that along its opposite direction are respectively given by 



p+{e)= I n(0)Cos^ ^" ^'^' d<i) (1) 



2{o-p) 



P-i9)= f U{(f>)Sin''^^—^d(j) (2) 

Jo 2 

where P-|_ (9) + P_ (9) — 1, and n' is any positive number. 

It is these probabilities P+ (9) and P_ (9) which constitute the actual observable quantities in our scheme which are 
determined by the distribution of spins 11 {(/)) of the particles emerging from the SR. The theoritical estimations of 
these probabilities crucially depend on how one calculates the quantity 11 (</)) whose evaluation, in turn, is contingent 
on the approach adopted for calculating the relevant time distribution 11 (t) mentioned earlier. Now, the important 
point to stress here is that the specification of such a time distribution is not unique in quantum mechanics. For 
the setup indicated in Fig.l, 11 (t) represents the arrival time distribution at the exit point (x =^ d) of the SR, which 
is also the distribution of transit times[t) as well as the distribution of times of interaction with the magnetic field 
within the SR. Thus, n(i) determining 11(0) is the key quantity which needs to be evaluated first. In the following 
section we focus on discussing the Bohmian procedure for calculating Ii{t) leading to Ii{(t>)- 

IV. THE EVALUATION OF n (i) AND n(0) USING THE BOHMIAN MODEL 

In the Bohmian model[3, l5|, each individual particle is assumed to have a definite position, irrespective of any 
measurement. The pre-measured value of position is revealed when an actual measurement is done. Over an ensemble 
of particles having the samewave function tp , these ontological positions are distributed according to the probability 
density p — \ij}\'^ where the wave-function ijj evolves with time according to the Schrodinger equation. For the purpose 
of the present article, the relevant key postulates of the Bohmian model can be expressed as follows: 

i) An individual 'particle' embodying the innate attributes like mass, charge, and the magnetic moment has a 
definite location in space at any instant, irrespective of whether its position is measured or not. In particular, the pre- 
measurement value of position is taken to be the same as the value revealed by a measurement of position, whatever 



be the measurement procedure. However, this feature does not hold for any other dynamical variable, such as the 
momentum or spin. 

ii) IV'pdx is interpreted as the probability of a particle to be actually present within a region of space between x 
and x+dx, instead of Born's interpretation as the probability of finding the particle within the specified region. 

iii) Consistent with the Schroedinger equation and the equation of continuity being interpreted as corresponding 
to actual flow of particles with velocity v{x, t), the equation of motion of any individual particle is determined by the 
following equation 

J{x,t) 1 dS{x,t) 

Here J(x,t) — p(x,t)v{x,t) is interpreted as the probability current density of actual particle flowaiiA the wave 
function is considered in the polar form ij: — Re^^/^ with p = R^ where R and S are the real functions of position 
and time. In the context of our SR example, coming to the question whether there can be predictions obtained from 
the Bohmian model which may allow to empirically discriminate it from the standard scheme, we analyse this isue as 
follows. 

Here it should be noted that the Schroedinger expression for the probability current density is inherently nonunique. 
This is seen from the feature that if one adds any divergence-free or constant term to the expression for J{x, t), the 
new expression also satisfies the same equation of continuity derived from the Schroedinger equation. However, it 
has been discussed by Holland[6j, followed by others 
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32| . that the probability current density derived from the 



Dirac equation for any spin-1/2 particle is unique, and that this uniqueness is preserved in the non-relativistic limit 
while containing a spin-dependent term in addition to the Schroedinger expression for J{x, t). Accordingly, a modified 
Bohmian equation of motion for the spin-1/2 particle has been proposed [6] which is of the form given by 

v(x, t) = — {\7S{x, i) + V log p{x, t) X s(t)) (4) 

m 

where s(t) — §x(i)'5^X^(^); x{t) is the time evolving spin state and the components of a are the Pauli matrices. 
The modified Bohmian equation of motion given by Eq.Q reduces to the original Bohmian form given by Eq.(3) 
if [VS*! 3> § jVlogpl, i.e., when the modulus of the spin-dependent term is negligibly small compared to the spin- 
independent term. In our setup, the pertinent parameters (viz. the magnitude of the magnetic field, the initial width 
and the peak velocity of the wave packet) are taken to be such that this condition is satisfied. Hence, throughout this 
paper, we will neglect any effect of the spin-dependent term in the modified Bohmian equation of motion. 

Next, we recall the general question of equivalence between the standard version and the Bohmian model of quantum 
mechanics in predicting the observable results pertaining to any Herniitian operator that has been much discussed[2|. 



along with some controversy |48| . Nonetheless, the central point we may stress that, given any example, if and only if 
the procedure for calculating an observable quantity is unambiguously defined in both the standard and the Bohmian 
versions of quantum mechanics, the very formulation of the Bohmian model guarantees the equivalence(at-least, in 
the non-relativistic domain) when the (common) formalism is applied. 

Now, coming to our example, for the initial state given by '^ {x,t = Q) = tp {x,t = 0) (E)x{t = 0), the time evolution 
within SR is subject to the interaction Hamiltonian flb.B. Note that as the spins of neutrons interact with the 



time-independent constant magnetic field, in order to conserve the energy, there will be changes in the momenta 
of the spatial parts associated with both the up and down spin components - these changes occur according to 
the potential energy of the spin-magnetic field interaction that develops corresponding to the up and down spin 
components respectively. 

In the setup considered here, the relevant parameters(viz. the initial momentum, mass, magnetic moment, and the 
magnitude of the applied constant magnetic field) are chosen to be such that the magnitude of the potential energy 
(wec/icr.B) that arises because of spin-magnetic field interaction is exceedingly small{K, 10~^ times) compared to the 
initial kinetic energy of the neutrons. This amounts to ensuring the magnitude of the momenta changes within the 
SR to be negligibly small compared to the initial momentum. Consequently, the time evolutions of the spatial and 
spin parts can be treated independent of each other, with the spatial wave function evolving freely within the SR(for 
the relevant justification based on a rigorous treatment, see Appendix A ). 

With respect to the setup shown in Fig.l, the initial wave function is given by 



il}{x,t^Q) = - -YTjexp 



" ~ 7T "I IfCX 



(5) 
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where the wave number k = mu/h, u being the group velocity and CTq is the initial width of the wave packet. Now, 
since the spatial wave function is considered to evolve freely within the SR, the freely propagating wave function at 
any instant t is given by 



^^"'')^(2^IfjVJ*^"P 
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+ ik{x~ ut/2) 



(6) 



where At = (Tq[\ 
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Writing Eq.® in the polar form tp{x,t) = _R(a;, t)e*'^(^'*)/'', and using Eq.(3), the Bohmian trajectory equation in 
this setup for an i-th individual particle having an initial position 2:^(0) is given by 



X^ (t) =ut + X^ (0) Vl + Pf^ (7) 

with (3 = A 2 i i and the subscript 'i' is used as the particle label. 

The Bohmian velocity of an z-th individual particle can then be written as a function of the initial position by using 

Eq.Q. The resulting expression is given by 

dx, jt) xoil3t 

vdxo.,t) = ^^^u+-^== (8) 

From the above equation it follows that only those particles initially located around the peak position of the Gaussian 
wave packet at xot ~ follow the Newtonian equation of motion. On the other hand, the particles initially located 
within the front-half( xoi — +ve) are all accelerated, while the particles initially lying within the rear-half( xoi ~ ~ve) 
are all decelerated [49]. It can then be seen that there will be be turning points of the Bohmian trajectories provided 
the following relation is satisfied between the instant t' of the turning point and the Bohmian initial velocity u, given 

by X0^|3t' 

u = — = (9) 

where a;oi is essentially —ve. For the purpose of discussions in this article, we consider the relevant parameters [u 
and ao) to be such that the above condition is not satisfied, thereby ensuring that, in our example, the particles do 



not have any turning pointf5Cf|. Next, in order to proceed with the Bohmian calculation of 11(0), we focus on the 
following two quantities pertaining to neutrons which are initially localised, say, within an arbitrarily chosen region 
between Xi{0) and Xi{0) + dxi{0). 

First, we note that the probability of neutrons to be actually localised between Xi{0) and Xi{0) + dxi{0) is given by 
\'4'{xi{0),t = 0)\^dxi{0). Secondly, we consider the probability of these neutrons to cross the region of the magnetic 
field enclosed in the SR within the time-interval, say, T and T + dT; this is denoted by IlB{T)dT where IIb{T) is the 
arrival time distribution at the exit point (x ~ d) of the SR, T being the time required by the i-th neutron to reach 
the exit point (x — d) of the SR. 

Given the notion of trajectories in the Bohmian model, and taking that all the neutrons initially localised within 
Xi{Q) and Xi{Q)+dxi{Q) ultimately pass through the entire region of the SR (i.e., all of them cross the exit point(x — d) 
of the SR, and none of them reenters) , the above two quantities can be equated so that 

TlBiT)dT - \^{x^{0),t^ 0)\^dx^{0) (10) 

whence 

UB{T)^\i;{xm,t^O)\'^ (11) 

Note that Eq.(ll) is a very general relation that determines the arrival time distribution in any causal trajectory 
model, valid for any form of the initial wave packet. In the particular case of a freely evolving Gaussian wave packet 
we are considering here in the context of the Bohmian model, we proceed as follows: 

By using Eq. ^, and replacing t by T, the quantity ^'j, ' on the right hand side of Eq.(ll) is calculated by writing 
Xi{0) in terms of Xi(T). Further, by invoking the earlier outlined interpretive ingredients (6) and (d) of the Bohmian 
model, we write the quantity \t/;{xi{0),t = 0)^ in terms of Xi{T). This is done by using Eq.(3) for ijj{x,t — 0), 
replacing the argument 'x' of this function by 'a;i(0)' for varying 'i' corresponding to different initial positions of the 
neutrons. Then, putting Xi(T) — d, using both Eq.Q and the Schroedinger expression for the probability current 
density J{x,t), it can be verified that Eq. lfTl]) reduces to 

UBiT) = J{d,T) (12) 

where 

{d-uTf\ r {d-uT)h^T 



is the current density of neutrons at the exit point (a; = d) of the SR. Here ax is the width of the wave packet at the 
instant T. Note that J{d,t) given by Eq. P^ is essentially a +ve quantity. Thus, Eg. p^ shows that the Bohmian 
arrival time distribution at x — d is the same as that given by the probability current density(henceforth designated 
as PCD) approach[23; i.e., 

IlBiT)=IlpcDiT) = Jid,T) (14) 

Then, comes an important element in our example, viz. the consideration of the time evolution of the spin state 
according to an equation that is obtained by decoupling the position and the spin parts of the Pauli equation (adopting 
the 'approximation' stated earlier whose justification is indicated in Appendix A), given by 
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IM7.B X = ih-^ (15) 

Now, we recall that in this particular example, the initial spin of any neutron is taken to be polarized along the 
+x-axis; i.e., 

x(0) = K). = ^[|t>. + li)J (16) 

Then, since the constant magnetic field B = Bz within the SR, confined between x — and x — d, is switched on at 
the instant (i = 0) the peak of the wave packet reaches the entry point (x — 0)of the SR, the time evolved neutron 
spin state x (t) under the interaction Hamiltonian H = fla.B is given by 

X (r) = exp (^^) X (0) = ^e^ ^ [|t), + e^^"^ |i) J (17) 

where uj — fiB/h and r is the spin- magnetic field interaction time; i.e., the time during which the neutron spin state 
evolves according to the spin-magnetic field interaction term. Eq.([T7]) implies that after a time interval r, the spin 
orientation of a neutron is rotated by an angle with respect to its initial spin-polarised direction where 



2lut (18) 



Subsequently, one may put r = T in Eq.(|18p. and use Ea. p4|) to obtain the probability distribution of spin 
orientations (n((/))) of the neutrons emerging from the SR. Then, as applied to our example, the equivalence between 
the Bohmian scheme and the probability current density approach would seem to follow; i.e., 

nB(0) = npcD(0) = J(d,<^) (19) 

V. SUBTLETIES IN THE BOHMIAN CALCULATION OF n(i) AND n(0) 

A. Let us focus on a critical feature that is involved in obtaining the expression for IIb{T) given by Eq. (|12p within 
the Bohmian framework. This requires using the equality given by Eq. ljlOl) . taking it to be true for all values of the 
initial positions Xi{0) of the neutrons, including even for those neutrons initially located outside the SR. Consequently, 
Hb (T) given by Eq. (|12p is the arrival time distribution at a; = (i(the exit point of the SR) pertaining to the entire 
ensemble of neutrons, including the neutrons initially located outside the SR. 

Then, since writing the probability distribution 11(0) — J{d,(j)) from IIb{T) of Eq. ([T2| using Eq. (fT8l) is contingent 
upon taking T ~ t, this means regarding the arrival time distribution at a; = d to be essentially the distribution 
of interaction times r over which the spin states of the neutrons belonging to the entire ensemble evolve under the 
interaction Hamiltonian fla.H. 
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Now, if the above proviso is to be satisfied, this would imply that, within the Bohmian model, even if a neutron 
is located outside the region of the SR at the instant (t = 0) when the magnetic field within the SR is switched 
on, the spin state associated with that neutron would have to be regarded as evolving, under the spin-magnetic field 
interaction, from the instant t = itself. However, a fundamental feature of the Bohmian model is that an individual 
particle such as a neutron is regarded as a localised entity in an objective sense, that has a definite location in space 
at any instant, and which embodies the innate properties such as the mass and the magnetic moment. Further, 
a relevant crucial point is that the spin-magnetic field interaction term jib.'B has the localised particle attribute, 
magnetic moment(/i), as the coupling constant. It is, therefore, arguable that in applying the Bohmian scheme to 
our example, the time evolution of the spin state associated with an individual neutron is to be considered subject 
to the interaction /icr.B only if the neutron is inside the SR within which the magnetic field is confined. Thus, from 
this point of view, the earlier mentioned proviso would be inadmissible for the neutrons which are initially located 
outside the SR. In other words, within the Bohmian model, this would mean taking the ontological position variable 
to be more fundamental than the global wave function in determining the instant from which the individual particle 
is subjected to the relevant interaction which is mediated through a localized property of the particle (in this case, 
the magnetic moment). 

Consequently, it would not be legitimate to take the Bohmian arrival time distribution (IIb{T)) at the exit point 
(x = d) of the SR to be the distribution of interaction times r for all the neutrons belonging to the entire ensemble. 
This would then mean that the Bohmian probability distribution of spin orientations of neutrons emerging from the 
SR cannot be written as J{d, (p), thereby implying in the context of our setup, an alternative route for calculating the 
quantities HB{t) and WBi't')- 

Such an alternative Bohmian procedure for calculating \lB{t) and 'ilB{<t') would entail the viewpoint that the spin 
state of an individual neutron which is initially located outside the SR would evolve under the spin-magnetic field 
interaction that acts using the neutron's localised magnetic moment as the coupling constant only after that neutron 
reaches the entry point {x = 0) of the SR. Therefore, for such a neutron, its arrival time at the exit point(a; = d) of 
the SR would not be the time over which its spin-magnetic field interaction occurs in passing through the SR. On 
the other hand, for a neutron initially located inside the SR, its spin state would, of course, start to evolve from the 
instant t = itself when the magnetic field is switched on. 

Consequently, in this approach for calculating the Bohmian probability distribution 11^ (</'), the time over which 
the spin of a particle evolves under the magnetic field enclosed within the SR is essentially the particle's transit time 
within the SR [not the arrival time at x = d). which is the time taken by it to cross the region between x = Q and 
X — d. Thus, in such a procedure, it is the computation of the transit time distribution within the SR (in contrast to 
the arrival time distribution at x ~ d computed in the earlier discussed Bohmian procedure) that would provide the 
prediction of the Bohmian distribution IIb{4>)- Prima facie, it is thus an open question whether the quantitiy n^((/)) 
calculated in this alternative way would agree with the result obtained(Eq.(23)) from the Bohmian calculational 
scheme discussed in the preceding section - a procedure which, in a sense, ascribes more fundamental ontological 
status to global wave function associated with an individual particle than to its ontological position variable. 
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B. An important point to be stressed here is that the calculational scheme discussed in Section V is specifically 
couched in terms of a Gaussian wave packet. In particular, the crucial expression given by Eq.(12) is obtained from 
the generic Eqs. (10) and (11) in a way that involves an explicit use of the Gaussian wave packet. Therefore, it 
is again another open question whether the equivalence between the Bohmian and the probability current density 
based approach in calculating the quantity Hsit) would hold good if one considers non-Gaussian wave packets. For 
example, one may consider a non-Gaussian wave function in the momentum space given by 



Hk) 



N 



{2 7ral) 
where k is the wave number and iV = ( 1 



■ exp 



(fc - fco)' 



4-?. 



1 + a sin 






(20) 



^(1 — e~"5")j is the normalisation constant. The above function 
contains four real parameters CTfc,fco, a and /3 among which ak is positive. The salient feature of the corresponding 
wave-packet, i.e. |(/)(fc)p, is its infinite tail with the probability of finding particles negligibly small outside a bounded 
region determined by the parameter cr^, while the wave packet is asymmetric and it reduces to the Gaussian form 
upon continuous decrement of |a| to zero. Note that the infinite tail is in contradistinction to other non-Gaussian 
forms found in the literature which are generated by truncating the Gaussian distribution. The asymmetry of the 
wave packet due to the sine function entails a difference between the mean and peak values of the wave packet, and 
deprives fco from appropriating either of these values. Similarly, the parameter ak does not denote the standard 
deviation of |(/)(fc)p. 

Using Ea. (j20p . the initial wave function in position space is just the Fourier transform of 4>{k) and is given by 



ip{x,t^O) 



1 



/2 7r 



(i){k)e'^'' dk 



N 



.1/4 



exp 



4 a' 



+ ik^x 



1 + iae 16 sinh 



n X 



4:0- 



(21) 



(2^^2)1 

where a = {2ak)~^. 

It would thus be an interesting exercise to study the arrival/transit time distribution pertaining to the freely evolving 
wave packet corresponding to Eq. pop and use the results obtained to probe the equivalence between the Bohmian 
and the probability current density based scheme in calculating the empirically testable probabilities relevant to our 
setup. 



VI. SUMMARY AND OUTLOOK 



The setup analysed in this article seems to provide a particularly interesting example in the context of the founda- 
tions of quantum mechanics since it is capable of empirically discriminating between not only the various suggested 
standard quantum mechanical schemes for calculating the arrival/transit time distribution, but it may also enable 
to empirically verify or falsify the specific unique predictions obtained from any realist trajectory model of quantum 
mechanics like that of the Bohmian type. This is irrespective of the question as regards which standard quantum 
mechanical scheme is the empirically correct one for calculating the observable probabilities in our example. While 
comprehensive computations are being pursued along the lines A and B indicated in Section V, here it may be useful 
to summarise the implications of the different possible outcomes of the relevant experimental study of such a setup: 
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a) If the experimental results turn out to corroborate the prediction obtained from IIb{4>){— ^pcd{4')) given 
by Eg. pop , this would not only validate a particular quantum approach for calculating the arrival/transit time 
distribution with its justification provided by the Bohmian model, but it would also mean that if one analyses our 
example within the Bohmian framework, one would have to infer that although the magnetic field is confined within 
the SR, it would have to affect, from t = 0, the time evolution of the spin state of even those individual neutrons that 
are initially located outside the SR.Such an effect within the Bohmian model is rather curious essentially because of 
the assumed objective (i.e., independent of measurement)localisation of an individual particle, and that the neutron's 
spin-magnetic field interaction is mediated through the localised particle attribute, viz. magnetic moment, as the 
coupling constant. Note that this type of question does not arise within the standard version of quantum mechanics 
because of the absence of the notion of objective localisation of an individual particle. 

b) On the other hand, if the alternative Bohmian calculational procedure indicated in A of Section V yields the 
probability distribution Il'g{(p) that is different from IIb{4>), and if the experimental results are in conformity with 
n'^(0), then, within the standard framework of quantum mechanics, this would pose a challenge to find out which 
particular quantum scheme, in this case, would yield results in agreement with the Bohmian prediction 11^ ((/>). 

c) Another possibility, apart from that the experimental results may not agree with either IIb{4') or 11^(0) for the 
Gaussian wave packet, stems from the line of study indicated in B of Section V - i.e., if by using non-Gaussian wave 
packet, it is found that the Bohmian prediction in our example disagrees with that obtained from the probability 
current density based quantum approach, this would again require studying which particular quantum mechanical 
approach for calculating the arrival/transit time distribution, as applied to our example, would lead to results in 
agreement with the relevant Bohmian prediction using non-Gaussian wave packet. 

Thus, whatever be the experimental outcome, nontrivial questions would arise. Further, in view of a number 
of insightful variants [7H10| of the Bohmian model that have been proposed using the notion of objectively defined 
trajectories of individual particles, it should be instructive to analyse the example formulated in this paper in terms 
of such models in order to find out the extent to which an empirical discrimination is possible. Such studies would 
reinforce the importance of pursuing a comprehensive experimental programme along the direction that has been 
indicated here. 

Appendix: A 

The Gaussian wave packet considered in the paper can be regarded as made up of plane wave components. The 
treatment given here pertains to an individual plane wave component. 

Let the initial total wave function of a particle be represented by 'ifi = ipQ (E) x, where -00 = Ae^^^ is the spatial 
part which is taken to be a plane wave with wave number fc, and x = {—7s{x+z + X-z) is the spin state polarized in 
the +x direction where x+z and x-z are the eigenstates of ct^. Now, we consider that the particle passes through a 
bounded region(called SR) that contains constant magnetic field directed along the -|-z-axis. 

The interaction Hamiltonian is Hint — /iO'.B where /i is the magnetic moment of the neutron, B is the enclosed 
magnetic field and a is the Pauli spin vector. Then the time evolved total wave function at i = r after the interaction 
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of spins with the uniform magnetic field is given by 

iHt 
*(x,r) = exp( ^)*(x,0) 

= -j^[ij+{^,T)®X+z+i^-{^,T)®X-z\ (A.l) 

I ^^A ■ ■ 

where "0+ (x, r) and ■(/;_ (x, r) are the two components of the spinor i/' = which satisfies the Pauh equation. 

We take the enclosed constant magnetic field as B = Bz. The time evolutions of the position dependent amplitudes(V'+ 
and V'- ) of the spin components of the wave function are, thus given by 



Eqs.(A.2) and (A. 3) imply that when a neutron having spin up interacts with the constant magnetic field within 
the SR, the associated spatial wave function (V'+) evolves under a potential harrier that has been generated due to 
the spin- magnetic field interaction, while for a neutron having spin down, the associated spatial wave function (■(/'-) 
evolves under a potential well. 

Here we will specifically consider the situation where E > |/ii?|. This is because, even if one uses low energy or 
ultra-cold neutrons having kinetic energy of the order of 5 x 10~^eT^, if the potential energy term (|/i-B|) has to exceed 
the kinetic energy term, one will require the magnitude of the magnetic field to be exceedingly high, to be of the 
order of lOT. Magnetic fields of such high intensity are difhcult to produce in the usual laboratory conditions, and 
therefore for all practical purposes, it suffices to treat the case where E > |/i-B|. 

Now, since V'- evolves under a potential well confined between a; = and x ^ d the reflected and transmitted parts 
are respectively given by 

/ — A —ikx V 1/V ^ ) ^ \ A\ 

"^^ " (fc + fcl)2 - (/fc - /fci)2e2»fcid ^ ' 



7/;- - 4pifc:r^ /A r\ 

^^ ~ {k + fci)2 - (fc - kire^ik^d ^A-^i 



where k — ^^^-E ^ j^-^ = v '"'■^ — tLJ. and d is the width of the SR arrangement which contains the uniform magnetic 
field. 

On the other hand, for ■0+ which evolves under a potential barrier, the expressions for the transmitted and the 
reflected part are written by replacing all the /ci's in Eas. (jA.4p and (JA.5P by k2 where k2 — '"^ — ^^—^ 

^^ " ^ (fc + k2f - {k - k^Ye^^^-d ^^■'^) 
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Ak-lcr,p-i'"i ik2d 

^^~ (fc + A;2)2-(fc-fc2)2e2''=^'^ ^ ' 

Now, we note that the solutions in these cases, given our initial state, consist of a reflected part traveling in the 
— Jc-direction and a transmitted part of the wave function, that travels in the +2;-direction. However, the reflected 
part of the wave function exists only to the left of the SR, and the transmitted part exists only to the right of the 
SR. Our objective here is to calculate the observable distribution of spins of the particles emerging from the SR. For 
this, we need to look at only the transmitted part of the wave function. Therefore, the final time evolved transmitted 
state that is relevant to our purpose is of the following form that represents an entangled state involving the spin and 
the spatial degrees of freedom 

*/ = -t^O^tX+z + ^tX-z) (A.8) 

where N is the normalized constant that can be written as N = f (V't^P + IV'tP)'^^- 

It is then seen that, in the regime E > \fJ,B\, using Eqs.(A.5)and (A. 7), one can rewrite Eq.(A.8) in the following 
form 

Apikx 

*/ = -j^{Ce"^'x+z + De'*^x~z) = i'oxW (A.9) 

whereby Eg. dA.Op is the form of the Larmor precession relation which is valid under the condition E > \fiB\ that 
has been calculated in terms of the explicit time evolved solutions of Eqs.(A.2) and (A. 3). 
Here 



C 



^ Re{xpj,)^ + Imiip:^)^ (A.IO) 



and 



D = ^Re{tlj+y + Im{%jj+Y (A.ll) 



^,=.an-^ (A.12) 



Note that, by using Eq. ljA.sp . one can find that 

8kki{P + kf)sin{kd)sin{kid) + 16k'^klcos{kd)cos{kid) 
''^^^' " (fc + A:i)4 + (fe - fci)4 - 2(fc + A:i)2(fc - kiycos{2kid) ^ ' ' 

8kki{k^ + kf)cos{kd)sin{kid) — 16k^kisin{kd)cos{kid) 
"^^^^^ " (fc + fci)4 + (fc - fci)4 - 2{k + fci)2(fc - ki)^cos{2kid) ^ ■ ' 

Similarly, by using Eq.(A.7), one can obtain the expressions for Re{ipj^) and Im{tp^) given by 

^ 8kk2{k'^ + k2)sin{ka)sin{k2d) + 16k^k2Cos{kd)cos{k2d) 

^^'^^' " (fc + ^2)4 + (fc - fc2)4 - 2(fc + fc2)2(fc - k2ycOs{2k2d) ^ ' ^' 
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, 8kk2{k'^ + k2)cos{ka)sin{k2d) — 16k'^k2sin{kd)cos{k2d) 

Im{il^T) = (fc + k^Y + (^ _ k2Y - 2(fc + k2Y{k - k2Ycos{2k2d) ^ -^^^ 

Now, let us examine in what limit one can obtain the standard expression for Larmor precession. 

Considering the more stringent hmiting condition E » |/i-B|, it can be seen that when the kinetic energy term of 
the Hamiltonian is much larger than the potential energy term, one can assume that the time evolution of the wave 
function occurs effectively due to a very shallow well and a very low harrier. This situation would correspond to the 
entire wave being transmitted while picking up just a phase. 

From the expressions for fc, fci, fc2 given earlier, in the limit E >> |/i-B|, we find that fc « fci « ^2- Then, in order 
to get the standard expression for Larmor precession, one can first set k — ki = k2 va Eqs. (jA.14p . (|A.15p . (|A.16p . and 
in (|A.17p . except when they appear inside the sine or cosine functions, since the latter terms are much more sensitive 
to the differences in the values of k, ki,k2- Eqs. (JA.14P and (|A.15P would then simplify to 

Reiip^) — sin{kd)sin{kid) + cos{kd)cos{kid) (A. 18) 



Im{'ipj,) — sin{kid)cos{kd) — cos{kid)sin{kd) (A. 19) 

Using the above expressions in Eqs. ()A.10|) and (jA.12p . we find that C = 1 and 0i — (fci — k)d. Similarly, rewriting 
equations (IA.16|1 and (|A.16p . and using Egs. ljA.Up and (|A.13p . we get D = 1 and 02 = (^2 - k)d. Therefore, Eq.(A.9) 
reduces to the form 

^f = —^ U^^^-^^'^X^y- + ^"'^^^'^^'^x-z) (A.20) 



Now. remembering that k — ^^™"^ , k\ = ^ f, — '~^ ^'^'^ ^2 — n — ^~~^' ^^ '^^^ binomially expand k\ and ^2 

around fc, and keep terms up to the order of ^, since we have already stipulated the condition k Ri ki ~ A;2. Then, 
{ki — k)d = —k^d, and using v — ^, we can write (fci — k)d — ~^f = wi. 
Similarly, (^2 - k)d = ^^ = ~^t. Therefore, we can write Ea. (|A.20p as 

Apikx 

^f = ^=- (e— *x+. + e^"*X-.) (A.21) 

where (j) — 2wi. The spin part of this equation is Eq.(17) given in the text, written by considering the time evolutions 
of the spatial and the spin parts to be independent of each other, with the spatial wave function evolving freely within 
the SR. 

It is important to stress here that the rigorous treatment given above of the time evolution of the quantum state 
of a neutral spin-1/2 particle in a constant magnetic field reveals that the standard expression for Larmor precession 
given by Eq.(14) in the text is essentially justified in the limit where the kinetic energy term is much higher than the 
potential energy term. Thus, subject to this condition being satisfied, the assumption of mutual independence of the 
evolutions of the spatial and the spin parts of the total wave function that has been used in the text is justified for 
the setup considered in this article. 
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